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One of the biggest challenges associated with a nulling interferometer-based 
approach to detecting extra-solar Earth-like planets comes from the extremely 
stringent requirements of pathlength, polarization and amplitude matching in 
the interferometer. To the extent that the light from multiple apertures are not 
matched in these properties, light will leak through the nuller and confuse the 
search for a planetary signal. Here we explore the possibility of using the co- 
herence properties of the starlight to separate contributions from the planet and 
nuller leakage. We find that straightforward modifications to the optical layout 
of a nulling interferometer will allow one to measure and correct for the leakage 
to a high degree of precision. This nulling calibration relaxes the field match- 
ing requirements substantially, and should consequently simplify the instrument 
design. 



Subject headings: techniques:interferometric 
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1. Introduction 

One suggested method for finding Earth-like planets orbiting other stars involves a space- 
based nulling interferometer operating in the thermal IR (Bracewell 1978; Woolf & Angel 
1998; Noecker 1999). The National Aeronautics and Space Agency has been considering 
such a mission: "Terrestrial Planet Finder" (TPF; Coulter, 2004), while the European Space 
Agency is planning an equivalent mission: Darwin (Kaltenegger et al. 2003). The primary 
challenge in detecting such planets is the extreme contrast ratio between the planet and its 
primary star (~ 10 6 ~ 7 at a wavelength of 10/im) and close angular separation (~ 0.1 arcsec). 
To prevent the planet light from being overwhelmed by photon noise from the primary, some 
way of selectively nulling the starlight is required. However, the extreme nulling ratio desired 
sets very stringent limits on the allowable instrument performance, specifically in terms of 
pathlength error (5<f> 2 ), polarization (5a 2 ) and amplitude mismatch (5a 2 ), all of which cause 
"leakage" in the null. Lay (2004) has shown that for certain nulling configurations, in 
addition to simple leakage terms of the form 5(f) 2 or 5a 2 , second-order coupling effects of 
the form 5(f)5a lead to stability requirements on the order of ~ 1.5 nm in path and ~ 0.1% 
in amplitude over hour-long timescales. These performance levels are almost an order of 
magnitude more stringent than previous estimates, and pose severe challenges to the TPF 
mission design. 

In this paper we explore a simple way of determining the nulling leakage based on 
exploiting the fact that the leakage light is coherent with the starlight, but not with the planet 
light. By making a simple modification to the interferometer back-end it becomes possible 
to measure the leakage terms and hence calibrate the null. This idea is an adaptation of the 
"Synchronous Interferometric Speckle Subtraction" concept proposed by Guyon (2004) for 
coronographic instruments. In the Section 2 we describe the original Bracewell nuller concept. 
We then introduce the nulling calibration technique and derive its expected performance, 
as well as show simple simulations of the application. In Section 3 we introduce the more 
complicated dual-Bracewell architecture actually being considered for TPF, as well as show 
how the calibration concept could be applied to it. We discuss the effects this might have 
on the TPF system design in Section 4. 

2. The Single Bracewell Nuller 

— * 

Consider Figure 1. Light is collected from two apertures separated by a baseline B and 
brought to a common point such that the wavefronts from the two arms are exactly n/2 out of 
phase with respect to each other going into the beam-splitter, for a total post-combination 
phase shift of n. Recall that a beam-splitter introduces an additional 7r/2 relative phase 



Fig. 1. — A simple beam combiner used in the Bracewell concept 



shift between reflected and transmitted beams. For a monochromatic interferometer the 
initial n/2 is simply a path adjustment. However, in the more realistic case of a broadband 
interferometer a more sophisticated approach becomes necessary, typically in the form of 
"phase plates", i.e. one beam is passed through a set of carefully controlled thicknesses of 
glass arranged such that the wavelength-dependent indices of refraction combine to give a 
7r/2 phase-shift across a 10-20 % bandpass (Wallace et al. 2004). At the 50/50 beam-splitter 
the two beams are combined, and the intensities can be measured at the two complementary 
outputs. Note that the planet-star separation is unresolved by the individual apertures, and 
hence in this "pupil plane" interferometer, detection of the planet signal requires that it be 
modulated in some fashion; in the Bracewell configuration this is accomplished by rotating 
the interferometer about the axis pointing to the star. 

If we neglect time- dependence and polarization of the electric fields, we can write 

E l = A ai e i{<pl) (1) 
E 2 = Aa 2 e t(,l>2+n/2) (2) 

where A is the amplitude of the field seen by an ideal aperture, <f>i, (fi 2 , cli and a 2 represent 
small phase and amplitude mismatches between the arms of the interferometer (0 <C 1 and 
a ~ 1). The two output beams become 

E l + E 2 e in ^ 2 

* = ^f— (3) 

E d = ^Lt* (4) 
and we recover the usual interferometric fringe as 



h = \Eb 
A 2 



2 



= ~2 (°1 + fl 2 + 2aia2 COS (01 _ 02)) (5) 
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h = \Ed\ 
A 2 

= ~2 (°i + a l ~ 2a i°2 c os(0i - 4>2)) (6) 

The phase difference 0i — 02 has two terms: instrumental phase differences (5<p) and a 
geometric term. In the ideal planet-search situation 5<p is small, the instrument baseline is 
perpendicular to the direction of the star, and the planet is at a small angle such that light 
from it has an additional phase given by 

P = y B ■ s (7) 

where B is the baseline vector separating the apertures, and s is the star-planet separation 
vector, both projected on the sky. By rotating the interferometer baseline about the direction 
to the star, it is possible to produce a time-variable signal in the Id output coming from the 
planet, while keeping the star nulled. If the null were ideal (5a = a<i — a\ ~ and <50 ~ 0) 
there would be no stellar leakage to overwhelm the planet signal in the dark output Id- 
Unfortunately, I p /I s ~ 1CT 6 and hence we would require 5a, 5(f) < 1CT 3 . 

In the absence of near-perfect amplitude and phase control, a detector at the "dark" 
output of the nulling interferometer will see light from two sources: the planet, which will 
vary according to the rotation of the interferometer, and the leakage, which may vary on all 
timescales. To the extent that this variability occurs at frequencies coincident with the planet 
signal, it will cause confusion and loss of planet Signal-to- Noise Ratio (SNR). However, if it 
were possible to measure the leakage terms after the beam combination, it might be possible 
to significantly relax the control requirements, making the instrument more feasible. In such 
a situation it is no longer necessary to require that leakage levels be extremely stable; instead 
it is only necessary that that uncalibrated changes of the leakage (e.g. due to shot noise) be 
smaller than the planet signal. 

One might consider monitoring the "bright" output of the nuller lb, and use conservation 
of energy arguments to infer that any reduction of photon counts in lb implies a corresponding 
increase in Id- However, for the typical TPF observation, Jf, ~ I s ~ 10 6 phot s _1 and the 
photon noise will be too great (^/h » I p ) to determine Id to sufficient precision in a 
reasonable amount of time. This approach would also be extremely sensitive to any changes 
in background or detector gain, as variations at the part in 10 6 level would render the 
calibration useless. 

A different approach has also been suggested (Danchi et al. 2003), where one can take 
advantage of the fact that instrumental mismatches are typically in path-length rather than 
phase, and hence will have a known wavelength dependence; this should allow one to use 
wide-bandwidth data to solve for and subtract systematic errors. However, such multi-color 
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approaches are still vulnerable to e.g. polarization mismatch errors, and they will require 
that the nuller operate with very wide band-passes, placing stringent requirements on the 
achromatic design elements. 

2.1. Calibration of a Single Bracewell Nuller 




Fig. 2. — Calibration of a Bracewell null. Fractions of the bright and dark nuller outputs 
are recombined with a controlled relative phase shift (9); the intensities of the calibration 
outputs can be used to solve for amplitude and phase mismatches in the input beams. 

The dominant part of light that leaks through the null comes from the star. Therefore, 
if we mix (e.g. at a beam-splitter) the electric fields of the leakage with that of a separate 
reference beam, also from the star, fringes will form as long as the relative path delays are 
maintained to within the coherence length of the light (Figure 2). On the other hand, light 
from the planet is not coherent with the starlight (i.e. (E star E p i anet ) = 0, where the brackets 
indicate an average over time), and hence will not form fringes. If we take a portion of the 
two outputs (e.g. 50%) of the nulling interferometer and recombine them again with a fourth 
beam-splitter, the resulting electric fields of these "calibration" outputs are 



E b e^ + E d e^ 2 
Ed = 7f (8) 

E b e^/V + E d 
E c , 2 = ^ (9) 

where is any additional phase introduced into one of the arms of the calibration interfer- 
ometer. The corresponding intensities are 
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A 2 

I c ,i(9) = ( a l + a l + (°i - a \) cos(0i - 2 ) + 2aia 2 sin(0i - <j> 2 ) sin(0)) (10) 
A 2 

I c ,2{0) = {al + a\ - (al - al) cos(0i - 2 ) - 2aia 2 sin(0i - <f> 2 ) sin(0)) (11) 

A fringe pattern forms in the output of the calibration interferometer; this fringe pattern 
contains information about the amplitude and phase mismatches of the input beams. This 
information can be extracted in a straightforward manner, given that 6 is entirely internal to 
the instrument and under our control. If we measure the intensities I c twice, having added 
A6 = tt/2 to 6 for the second reads, and labeling the reads "A", "B" , "C" and "£>", i.e. 



A = 


lc,i{0 


= 0) 


B = 


h,2{e 


= 0) 


C = 


lc,l(0 


= 7T/2) 


V = 


ic,2(e 


= tt/2) 



(12) 
(13) 
(14) 
(15) 

these measurements can be used to "reconstruct" the output Id caused by amplitude and 
phase mismatches, i.e. the leakage through the null. 

h JA-BfHC-Vf 

This leakage can be subtracted from the measured I a, the remaining light is from the 
planet. While it is true that the planet light will also pass through the calibration interfer- 
ometer and create an interference fringe, because the light from the planet is not coherent 
with the starlight, yet both fall on the same detector, the fringe contrast of the fringe due 
to the planet will be reduced by a factor of I p /I s ~ 10 -6 and will have a negligible impact 
on the measurement of the leakage parameters. 

A nulling interferometer will have starlight leak through the null for two reasons: the first 
is due to amplitude, path (and polarization) mismatches as discussed above. The second is 
simply due to the finite size of the star. Even a perfect null only blocks light that is exactly 
on-axis; for typical stellar sizes and distances the star will subtend an apparent angle on 
the order of a milli-arcsecond, large enough that a non-negligible amount of light will leak 
through the null. This leakage cannot be removed by the proposed calibration technique. 
However, the level of this leakage is set by the length of the nulling baseline, and should not 
vary on timescales comparable to the planet signal. Nonetheless, shot noise from the light 
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that leaks through can overwhelm the planet signal; it is this leakage term that limits the 
maximum size of the nulling baselines. 

It is interesting to note that the reconstructed leakage through the null has more favor- 
able noise properties than a leakage term reconstructed from the the bright output alone; 
this is a manifestation of "heterodyne gain". Assuming that the calibration intensity mea- 
surements are dominated by photon noise (i.e. o\ = A) the uncertainty of the calibration 
measurement is given by 



E 

Ri=I b ,A..V 



A 2 (5a 2 + 



OR, 



fit. (17) 

(18) 
(19) 



In other words, the reconstructed leakage I d has the same Signal-to-Noise ratio as what one 
would find by measuring I d itself at the dark output, except that now there is negligible 
contamination due to light from the planet; the mixing process allows one to separate the 
two cleanly. 

For 5a ~ 10~ 3 , 5(f) ~ 1CT 3 , A 2 = tI s and I s ~ 10 6 phot s -1 , assuming r ~ 100 seconds 
one should expect to measure Id to a precision of ~ 4 x 10~ 3 . The price for this calibration 
is a loss in the effective throughput of the instrument, as half of the light from the planet 
gets lost in the calibration system. Clearly there is a trade-off between the need to minimize 
photon noise from the planet vs. noise due to time-variable null leakage; we will discuss the 
effect this has on the overall instrument design in Section 4. 



2.2. Simulations of a Calibrated Single Bracewell Nuller 

We have simulated the operation of a simple two-aperture Bracewell interferometer 
observing an Earth-like planet. We assumed a baseline of 50m, 5.5-m diameter collecting 
apertures, and an Earth-sized planet in a 1 A.U. orbit with an albedo of 0.3. The central star 
was modeled as a 5700K, 1 R & blackbody located at a distance of 10 pc. Total integration 
time is 2.16 x 10 5 seconds, spread out over 2 full 360-degree rotations of the array, with 
300 seconds of integration time at each orientation. Path-length errors for each integration 
were modeled as 10-nm rms Gaussian random noise, and amplitude mismatch errors had 
a fractional amplitude of 0.001. Photon noise was applied to all measured intensities. The 
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system efficiency was assumed to be 30%. We assumed the system to have five l//m wide 
spectral channels ranging from 8 to 12/im in central wavelength. 

It is clear from Fig. 3 that the leakage due to amplitude and phase mismatches com- 
pletely overwhelms the planet signal, while the calibration successfully removed the leakage 
and recovered the planet. However, the required integration time and aperture size was 
considerable - a consequence of the single-Bracewell design that is limited by the trade-off 
between null depth and angular resolution. With only two apertures, the finite size of the 
stellar disk will cause light to leak through the null if the baseline is long enough to resolve 
the planet. Nevertheless, we reconstruct images of the target system using the standard 
cross-correlation analysis described by Lay (2005). The simulated lightcurves (both raw 
and calibrated) are cross-correlated with a series of templates. The templates are a func- 
tion of the position of a planet, and are given by the theoretical response of a two-aperture 
interferometer, i.e. Eqn. 5. The results are shown in Fig. 4. 



TPF has been envisioned not as a simple two-aperture Bracewell design, but rather a 
somewhat more sophisticated dual-Bracewell design. This is for two reasons: the response 
function of the two-aperture Bracewell is symmetric, leading to a 180-degree ambiguity in 
the position angle of any planet (Figure 4), and in the two-aperture Bracewell the planet 
signal can only be modulated on the array-rotation timescale. That timescale is ~ 1000 
seconds, long enough to be susceptible to systematic errors due to instrumental drifts. 

A solution to both of these problems is achieved in the form of the "dual chopped 
Bracewell" design outlined in Figure 5 . In this configuration two Bracewell pairs are placed 
next to each other, and the nulled outputs from each are recombined with relative phase shifts 
that are switched between values and ir. Such phase chopping produces an asymmetric 
response on the sky, and since the internal phase shift can be adjusted at high frequency, 
it allows one to remove many forms of systematic error associated with the instrument gain 
and background. 

We derive the output of a dual Bracewell system below 



3. 



The Dual Bracewell Nuller 



Ei 
E2 
E 3 
E 4 



Aa 2 e l ^ /2) 

Aa 3 e^ s) 

Aa A e l ^ + ^ 



(20) 
(21) 
(22) 
(23) 
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(24) 

The apertures are combined pairwise to create "Left" and "Right" combinations as follows 

E 1 + E 2 e^ 

~ 7f ^ 

F El6tn/2 + E " (26) 

- — 71 — ^ 

F E " + E ^' 2 (27) 
" 71 (27) 

E 3 e-/ 2 + E 4 

E b , R - -j= (28) 

Combining the two nulled outputs with relative phases 0, tt, which we label the "A" and "B" 
states respectively, yields 

E dL + E d>R e lw/2 
Ami, A = j= 1/9) 



72 

E ddl , B = E ^ + E J_^ 12 (30) 



and the corresponding intensities are 



A 2 

hdi,A = [a\ + a\ + a\ + a\ (31) 

— 2a 4 a 2 cos(0i — 2 ) + 2aia 3 sin(0i — 3 ) — 2aia 4 sin(0i — 4 ) 

-2a 2 a 3 sin(0 2 - 03 ) + 2a 2 a 4 sin(0 2 - 4 ) - 2a 3 a 4 cos(0 3 - 4 )) 
A 2 

hdi,B = -j- (a? + a 2 + «3 + a 4 ( 32 ) 
— 2a!a 2 cos(0i — 2 ) — 2aia 3 sin(0i — 3 ) + 2a!a 4 sin(0i — 4 ) 
+2a 2 a 3 sin(0 2 - 3 ) - 2a 2 a 4 sin(0 2 - 4 ) - 2a 3 a 4 cos(0 3 - 4 )) 

The final quantities used in the image reconstruction are the "sine" and "cosine" chops. The 
sine chop is given by 

Isin = Iddl,A — Iddl,B (33) 

= A 2 (aids sin(0i - 3 ) - a 2 a 3 sin(0 2 - 3 ) - aia 4 sin(0i - 4 ) + a 2 a 4 sin(0 2 - 4 )) 

This gives a response that is asymmetric with respect to the phase center of the instrument, 
and hence will allow one to determine the position angle of any planet without the 180- 
degree ambiguity of the single-Bracewell configuration. However, note that in this case, small 



-10- 



amplitude and phase mismatches can couple as the sine chop now contains error terms of the 
form SaS(j). A particularly challenging aspect of this coupling is that there are higher-order 
terms that couple non-linearly. Given that these mismatches occur on a range of timescales, 
the non-linear mixing can inject noise at frequencies where the planet signal is maximized, 
even in the presence of servo-control systems that try to minimize the mis-matches. 

The "cosine chop" is 

Icos = Iddl,A + Iddl,B (34) 

A 2 

— — [a\ + a\ + a\ + a\ - 2a 1 a 2 cos(0i - <f> 2 ) - 2a 3 a 4 cos(0 3 - 4 )) 

In general the sine chop is more useful to the image reconstruction, as the cosine chop 
will include contributions from symmetric sources such as zodiacal light. In any case, if 
the instrument was switched between states A and B quickly enough to remove drifts in 
thermal background, gain, etc., then these quantities are primarily determined by the planet. 
However, there are still the effects of variable leakage to consider and avoid. 



3.1. Calibration of a Dual Bracewell Nuller 

The concept of coherent calibration using the non-nulled outputs from the nulling inter- 
ferometers can be applied in a straightforward manner to the dual Bracewell configuration 
(Fig. 6). In this configuration, the bright (non- nulled) outputs from the pairwise millers are 
combined, and this beam is then mixed with the combined, nulled output. The resulting 
electric fields and intensities are easily found 

E cl = dd ' 2 + bb ' 2 L^t (35) 

v2 

E dd , 2 e™/ 2 + E bb ^ e /V2 
E c2 = -j, (36) 

where is an additional, controllable phase shift that can be introduced between the bright 
and dark beams. The associated intensities are, after some tedious algebra, found as 

A 2 

hi,A = — (3 (a\ + a\ + a\ + a 2 ) (37) 

+2v / 2cos(6») {-a\ + a\ - a\ + a 2 ) 

+4v / 2cos(^) (aia 3 sin(0i - 3 ) - a 2 a 4 sin(0 2 - 04)) 

+4:^/2 sin(^) (— a\a 2 sin(0! — 2 ) — a i a 4 cos(0i — 4 ) 
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+a 2 a 3 cos(0 2 - h) ~ a 3 a 4 sin(0 3 - 4 )) 

—2a\a 2 cos(0i — 02) — 60103 sin(0i — 3 ) + 2aia 4 sin(0i — 4 ) 

+2a 2 a 3 sin(0 2 - 03) - 6a 2 a 4 sin(0 2 - 4>i) ~ 2a 3 a 4 cos(0 3 - 4 )) 
A 2 

I C 2,A = —^{al + al + al + al) (38) 

+2v / 2cos(6') {a\ -a 2 2 + a 2 3 - a 2 ) 
+4a/2 cos(^) (— aia 3 sin(0i - 3 ) + a 2 a 4 sin(0 2 - 04)) 
+4v^2 sin(0) (a±a2 sin(0i — <f> 2 ) + aia 4 cos(0i — 4 ) 
-a 2 a 3 cos(0 2 - 03) + a 3 a 4 sin(0 3 - 4 )) 

— 2aia 2 cos(0i — 2 ) — 6aia 3 sin(0i — 3 ) + 2aia 4 sin(0! — 4 ) 

+2a 2 a 3 sin(0 2 - 3 ) - 6a 2 a 4 sin(0 2 - 4 ) - 2a 3 a 4 cos(0 3 - 4 )) 
A 2 

I*,b = 777 (3 {a\ + a\ + a\ + a*) (39) 
lb 

+2^008(61) (- a \ + a 2 2 -a\ + a 2 ) 
+4v / 2cos(6 l ) (— aia 3 sin(0! — 3 ) + a 2 a4sin(0 2 — 4 )) 
+4v / 2sin(^) ( — a\d2 sin(0i — 02) + ai<2 4 cos(0i — 4 ) 
-a 2 a 3 cos(0 2 - 03 ) - a 3 a 4 sin(0 3 - 4 )) 

— 2aia 2 cos(0i — 02) + 6aia 3 sin(0i — 3 ) — 2aia 4 sin(0i — 4 ) 

-2a 2 a 3 sin(0 2 - 3 ) + 6a 2 a 4 sin(0 2 - 4 ) - 2a 3 a 4 cos(0 3 - 4 )) 
A 2 

I C 2,b = — (3 (al + a 2 2 + a 2 3 + a 2 4 ) (40) 
lb 

+2-\/2cos(0) (al -a 2 2 + a 2 3 - a 2 ) 
+4V2 cos(6) (a x a 3 sin(0i — 3 ) — a 2 a 4 sin(0 2 — 4 )) 
+4v / 2sin(6 l ) (aia2sin(0i — 2 ) — aia 4 cos(0i — 4 ) 
+a 2 a 3 cos(0 2 - 3 ) + a 3 a 4 sin(0 3 - 4 )) 

— 2aia 2 cos(0i — 2 ) + baia 3 sin(0! — 3 ) — 2ai<2 4 sin(0! — 4 ) 
-2a 2 a 3 sin(0 2 - 3 ) + ba 2 a 4 sin(0 2 - 4 ) - 2a 3 a 4 cos(0 3 - 4 )) 



If we measure the outputs I c i,I C 2 with 6 = and 8 = n/2 , in the fashion analogous 
to what is done for the single Bracewell case (i.e Eqns. 12 & lb), we can reconstruct the 
output from J^i; this works for both chop states. Hence we can recover I sin and I cos . 

The layout shown in Figure b is not the only one possible. In fact, one could use any 
bright output as an input to the calibrator, including the outputs from just a single nuller 
(E b>L or E b R ). As will be discussed in Section 5.4 there are some advantages to this particular 



- 12 - 



layout, in that it can be made fully symmetric. 

3.2. Simulations of a Calibrated Dual Bracewell Nuller 

We have simulated the operation of a standard 4-aperture dual-Bracewell interferometer 
observing an Earth-like planet (Fig. 8). Parameters were identical to the single-Bracewell 
case where applicable; the array was modeled as a linear 4-element array with a 50-m baseline. 
The integration time was split equally between the chopped states A and B. As in the single- 
Bracewell case, images were reconstructed using cross-correlations of expected signals as a 
function of planet position (Eqn 33 & Eqn 34). 

The gain from application of the coherent calibration approach becomes clear in Fig. 
9. Here we create model images for three cases of noise: 1, 10 and 50 nm of phase and 
0.01%, 0.1% and 0.5% of amplitude mismatch respectively. The calibration easily recovers 
the planet in all cases, whereas for the raw data only the lowest-noise case yields a planet 
detection in the image. 

4. Implications for TPF Instrument Design 

Using coherent calibration it is possible to measure starlight leakage through the null 
due to systematic path and amplitude mismatches to high precision. In the case where 
such leakage is the dominant noise source, this allows for a considerable relaxation in the 
associated path-length and amplitude control requirements. However, there are many sources 
of noise in the TPF instrument, and not all of them are amenable to this type of calibration 
(Fig. 10). For instance, local zodiacal light produces large fluxes on the detector at longer 
wavelengths; Poisson fluctuations in that flux can wash out the signal from a planet. 

Given that any TPF design must carefully weigh the relative contributions for many 
noise sources, the ability to calibrate and remove a large set of noise terms will undoubtedly 
influence what the optimal system design will be. A full investigation of such a system 
optimization is beyond the scope of this paper. However, given that the effects of small 
errors in phase and amplitude mismatch are significantly worse at shorter wavelengths, e.g. 
oc A~ 3 (Lay 2005) it is likely that a TPF-like instrument equipped with a nulling calibration 
system will be able to work at shorter wavelengths, something that may be particularly 
valuable given the presence of spectral features associated with H 2 around 6.3 /zm and 
CH 4 at 7.7/xm (Des Marais et al. 2002). Unfortunately, the underlying ~ 300K planetary 
blackbody emission is dropping rapidly shortward of 10 /iin, making observations in that 
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wavelength range very difficult, coherent calibration notwithstanding. 

Perhaps the clearest illustration of the effect of coherent calibration is apparent in 
Fig. 11, which shows the SNR values for calibrated and uncalibrated versions of a notional 
TPF configuration, as a function of the level of path-length stability achieved. Note that we 
considered only the amplitude-phase cross-term in this calculation; there are additional terms 
that a real instrument would have to consider (e.g. related to polarization). Clearly, if ~ lnm 
levels are possible, coherent calibration may not be worthwhile. However, this represents a 
fractional stability approaching 1 part in 10 4 , which may be impossible to achieve in practice, 
in which case coherent calibration is very useful. 

5. Limitations & Alternative Approaches 

There are a number of practical considerations and challenges to implementing this 
concept in a working interferometer. 

5.1. Low Fringe Visibility 

The leakage calibration quantities are extracted from a fringe in the calibration inter- 
ferometer. For a deep null (h/Id ~ 10 6 or more) this implies a very low fringe contrast or 
visibility (~ a/ hi Id)', such low fringe visibilities may be hard to measure in the presence of 
certain kinds of systemtic noise (e.g. detector noise) that can produce an additive bias to 
the measured fringe contrast. However, such additive biases can be measured and removed 
- and they are most easily measured with a brighter source. It is therefore always preferable 
to have the brightest possible reference beam (as long as the detector isn't saturated), even 
if the SNR of the leakage measurement is to first order independent of the brightness of the 
"bright" input to the calibration interferometer. 

5.2. Thermal Noise 

While thermal foreground emission (e.g. from the telescope optics) is usually considered 
incoherent in that it does not form fringes at the nuller, this is not the case for a calibrator. 
The calibration interferometer in effect forms a Mach-Zender interferometer in the instru- 
ment, and hence any light that enters one input port will form a fringe in the calibrator. 
We note that this fringe can be measured and corrected simply by blocking one input of the 
nuller and measuring the amplitude of a remaining fringe. It should also be noted that for 
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all but the longest wavelengths the amount of thermal emission that enters the system and 
forms a fringe in the calibrator will be substantially smaller than the flux from a planet. 

5.3. Splitting Ratio 

For the purposes of this paper, we have assumed that half of the "dark" output is split 
off and mixed in the calibration interferometer, with the associated loss of SNR. However, 
the optimal choice of splitting ratio is probably not 50/50 and may in fact be much smaller. 
Since the frequency at which the planet signal varies is on the order of the array rotation 
timescale, it may be sufficient to only split off enough light to measure changes in the nuller 
leakage on that much longer timescale. 

5.4. Non-singlemode Effects 

In this discussion we have treated the incoming electric fields as single-mode wavefronts. 
Such an approach is appropriate for initial explorations of the concept; it is also a reasonably 
good approximation in the case of a space-borne instrument equipped with single-mode filters 
of the type envisioned for TPF (Lay 2005). However, care may be necessary in designing 
the instrument so as to ensure that the calibration system samples the exact same spatial 
mode as the science detector. In practice there may be small alignment differences that 
cannot easily be removed. However, we suggest that these differences can be calibrated as 
follows. If the layout shown in Fig. 6 is made fully symmetric by adding a sampling beam- 
splitter in Eddi and mixing Ejdi with Ebbi, it then is possible to interchange the "bright" and 
"dark" outputs merely by adjusting the phase shifts applied to the inputs E 2 and £4 (which 
interchanges the bright and dark outputs on the pairwise millers.) The difference between 
the two resulting estimates of the leakage can provide a diagnostic that could be used for 
further calibration and/or system alignment. 

6. Conclusion 

We discuss the concept of coherent calibration of an interferometric null and its ap- 
plication to the Terrestrial Planet Finder instrument. We find that such an approach to 
calibration greatly relaxes the required levels of electric field matching, and associated sta- 
bility requirements. This should improve instrument performance, particularly at shorter 
wavelengths where the effects of path-length control limitations are most severe. 
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Fig. 3. — (Top) The raw intensities produced by the simulated two-aperture Bracewell 
interferometer, in a single wavelength channel, as a function of array rotation angle. (Bottom) 
The calibrated output, together with the ideal time series that would be expected from an 
interferometer with no noise. The ideal curve has been offset for clarity. Note that calibration 
has removed noise levels that were a factor of 10 3 times greater than the planet signal. 
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Fig. 4. — (Top Left) Reconstructed image using simulated uncalibrated multi-channel data 
from a single Bracewell nuller. (Top Right) Reconstructed image of the same system, now 
using the calibrated data. (Bottom Left) The image that would be expected from an ideal 
(noiseless) system. (Bottom Right) Reconstructed calibrated image, with the central psf 
feature (due to leakage through the null) divided out, using a rotationally- averaged raw image 
as the template. Note the 180-degree position ambiguity inherent in 2-aperture instruments. 




Fig. 5. — Schematic layout of the dual, chopped Bracewell configuration. The two aperture 
pairs 1,2 and 3,4 are combined so as to create nulls; the nulled outputs are then cross- 
combined. A relative phase of or n is imposed between combinations. 
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Fig. 6. — Schematic layout of the dual, chopped Bracewell configuration with calibration. 
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Fig. 7. — (Top) The intensities produced by the simulated four-aperture dual, chopped 
Bracewell interferometer, for a single spectral channel, as a function of array rotation angle. 
"Raw" is the intensity I S i n including simulated planet and noise sources as explained in the 
text. (Bottom) The calibrated output, together with the ideal time series that would be 
expected from an interferometer with no phase noise. The ideal curve has been offset for 
clarity. 
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Fig. 8. — (Top Left) Reconstructed image using simulated uncalibrated data from a dual 
Bracewell nuller. (Top Right) Reconstructed image of the same system, now using the 
calibrated data. (Bottom Left) The image that would be expected from an ideal (noiseless) 
system. (Bottom Right) The image reconstructed from the calibrated cosine-chop. The 180- 
degree position ambiguity disappears for the sine-chop data, as the dual Bracewell nuller 
provides the necessary phase information to determine the position of the planet. Also, note 
that the sine chop includes terms that depend on the phase difference between apertures 
1-3, 2-4, and 1-4, with longer effective baselines; hence the reconstructed image shows much 
higher spatial resolution. Note that there is only one planet in the simulation; the dirty map 
produced by the cross-correlation does include artifacts. A proper deconvolution, e.g. based 
on the CLEAN algorithm (Draper et al. 2006) is beyond the scope of this paper. 
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Fig. 9. — Simulated images reconstructed from raw and calibrated data, with three different 
noise levels. Top: f nm of phase noise and 0.01% amplitude mismatch. Middle: 10 nm 
of phase noise and 0.1% amplitude mismatch. Bottom: 50 nm of phase noise and 0.5% 
amplitude mismatch. 
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Fig. 10. — Selected noise terms in the TPFI signal-to-noise budget as a function of wave- 
length. Terms include local zodiacal emission, starlight leakage due to the finite size of the 
central star, and "systematic" leakage terms due to small path-length and amplitude mis- 
matches. We emphasize that these calculations are merely approximate; a full, high-fidelity 
SNR budget for a nulling interferometer is beyond the scope of this paper. Nonetheless, it 
should serve to illustrate the point that - as is appropriate in any well-optimized instrument 
system - there are multiple competing noise sources of similar magnitude. 
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Fig. 11. — The SNR of an Earth-like planet seen by a TPF Interferometer, for various level s 
of systematic noise, with and without coherent calibration. The calibration sacrifices photons 
in return for a much greater tolerance against systematic errors. 



